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ABSTRACT 

We  consider  the  problem  of  simultaneous  power  and  code-channel  allocation  of  a  secondary  transmit¬ 
ter/receiver  pair  coexisting  with  a  primary  code-division-multiplexed  (CDM)  system.  Our  objective  is  to 
find  the  optimum  transmitting  power  and  code  sequence  of  the  secondary  user  that  maximizes  the  signal- 
to-interference-plus-noise  ratio  (SINR)  at  the  output  of  the  maximum  SINR  linear  receiver  while,  at  the 
same  time,  the  SINR  of  all  primary  users  at  the  output  of  their  max-SINR  receiver  is  maintained  above  a 
certain  threshold.  This  is  an  NP-hard  non-convex  optimization  problem.  In  this  paper,  we  propose  a  novel 
feasible  suboptimum  solution  using  semidefinite  programming.  Simulation  studies  illustrate  the  theoretical 
developments. 


I.  Introduction 

Recent  experimental  studies  [1]  indicate  that  most  of  the  licensed  radio  spectrum  experiences  low  utilization. 
Cognitive  radio  (CR)  [2]  emerged  as  a  promising  technology  that  improves  spectrum  efficiency  and  utilization 
by  allowing  the  secondary  users/networks  to  share  the  spectrum  that  is  licensed  by  primary  users.  As  licensees, 
the  primary  users  have  always  higher  priority  to  use  the  spectrum  [3].  The  underlying  challenge  of  this 
technology  is  to  guarantee  the  Quality-of-Service  (QoS)  requirements  of  the  primary  system  yet  to  maximize 

QoS  for  the  secondary  users  [4] -[7].  Power  control  for  cognitive  CDM  systems  was  considered  in  [8]  where 
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the  proposed  method  does  not  take  full  advantage  of  signal  processing  at  the  receiver  end,  and  thus  restricts 
the  throughput  of  secondary  users.  Joint  beamforming  and  power  allocation  algorithms  for  cognitive  radio 
networks  were  presented  in  [9],  [10],  while  auction  mechanism  for  power  control  were  developed  in  [11], 

In  this  paper,  we  consider  a  secondary  system  with  a  code  division  multiplexed  (CDM)  mode  coexisting 
with  a  primary  CDM  system.  We  study  the  problem  of  designing  a  power  and  code-channel  allocation  scheme 
for  the  secondary  user  that  maximizes  the  output  SINR  of  the  maximum-SINR  linear  receiver  filter  under 
SINR  QoS  constraints  for  all  primary  users  and  a  peak  transmission  power  constraint  for  the  secondary  user. 
This  is  a  non-convex  NP-hard  problem.  In  this  paper,  we  propose  a  novel,  realizable  suboptimum  solution 
using  semidefinite  programming. 

The  rest  of  the  paper  is  organized  as  follows.  Section  II  is  devoted  to  CR  system  model  specifics  and  the 
formulation  of  the  optimization  problem.  In  Section  III,  we  present  our  proposed  solution  while  in  Section  IV 
we  propose  an  iterative  procedure  based  on  semifinite  programming  solution,  which  converges  to  a  suboptimal 
feasible  solution.  The  performance  of  the  proposed  scheme  is  evaluated  through  simulations  in  Section  V.  A 
few  concluding  remarks  arc  drawn  in  Section  VI. 

II.  System  model  and  problem  statement 

We  consider  a  primary  CDM  network  with  K  users  and  processing  gain  (code  sequence  length)  L.  We 
also  consider  a  secondary  system  (with  CDM  mode  of  operation)  in  the  uplink  spectrum  band  of  the  primary 
system  (Fig.  1).  For  simplicity  in  presentation,  we  assume  that  all  signals  propagate  over  a  plain  (no  multipath) 
additive  white  Gaussian  noise  channel.  We  denote  by  hi  and  qt,  i  =  1,2,...,  K ,  the  path  coefficients  from 
user  i  to  the  base  station  and  the  secondary  receiver,  respectively.  The  path  coefficients  from  the  secondary 
transmitter  to  the  base  station  and  to  the  secondary  receiver  are  denoted  by  hs  and  qs,  respectively.  All  path 
coefficients  are  modeled  as  Rayleigh  distributed  random  variables  that  are  independent  across  user  signals 
and  remain  constant  during  several  symbol  intervals  (quasi-static  fading).  After  chip-matched  filtering  and 
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sampling  at  the  chip  rate,  the  received  signal  at  the  primary  base  station  can  be  represented  as 

K 

r  =  ^  '  \/iw h jS i bt  T  \J Eshsssbs  T  r ip,  (1) 

i=  1 

while  the  secondary  received  signal  is 

K 

y  —  ''y  \  \J +  \J Esqsssbs  T  ns,  (2) 

i= 1 

where  Et.  bi  and  st  denote  the  bit  energy,  the  information  bit,  and  the  normalized  signature  vector  of  the 
primary  user  i,  respectively;  Es.  bs  and  ss  denote  the  bit  energy,  the  information  bit  and  the  normalized 
signature  vector  of  the  secondary  user,  respectively;  np  and  ns  represent  AWGN  Af(0.  a2l)  at  the  base 
station  and  at  the  secondary  receiver,  correspondingly. 

The  lineal-  filters  at  the  base  station  and  secondary  receiver,  which  exhibit  maximum  output  SINR,  can  be 
represented,  respectively,  as 

WMSiNR,i  =  cR^1Sj,  k  =  1, 2, . . . ,  K 
WMSINR,S  =  cR-1Ss, 

where  R  =  £'{rrT},  R  =  £'{yyT},  c,  c  >  0  and  E{ ■}  denotes  statistical  expectation.  The  output  SINR  of 
primary  user  i  at  the  base  station  and  the  secondary  user  at  the  receiver  can  be  expressed  as  in  (4)  and  (5) 
at  the  top  of  the  following  page,  where  R/j  and  R/s  are  the  “exclude  /”  or  “exclude  s ”  data  autocorrelation 
matrix,  defined  as 

K 

R/*  =  ^2  Ekhlsksl  +  Esh2sSssJ  +  cr2I 

k=l  ,k^i 

R/s  =  E  +  a2I. 

k= 1 

In  a  cognitive  radio  setup,  the  secondary  transmitter  has  to  guarantee  the  QoS  of  all  primary  users.  Hence, 
our  objective  is  to  find  the  transmission  bit  energy  Es  and  the  real-valued  signature  ss  that  maximize  SINRs 
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SINRi  = 


^'{lwMS/ATR,j(vAEi^i^si)|2} 


E{  I  wMS/iVR,i(Efc=i,fc^i  VEkhkSkbk  +  \/Ejhsssbs  +  np)|2} 


=  ^/ifsfR^s,:,  (4) 


SJJV.R,  = 


■^{lwMS/AT_R,s(V^s(/s&sss)|2} 

^{lwMS/iVR,s(Ef=i  VEkQkSkbk  +  ns)|2} 


=  Esq2s  sjR^Ss 


(5) 


under  the  constraints  that  SINRi,  i  =  1,2 , ...  ,K,  are  above  a  certain  threshold  a,  i.e. 

max  EljS^RT1  ss  (3) 

Es,Ss  s  /s 

subject  to  EihfsJUJ^  s i  >  a,  i  =  1,2,...,  K 

S^q  —  1  /  •  <f  /  ■ 

as  —  -*-5  -L/s  _  ^max 


where  Emax  denotes  the  maximum  bit  energy  for  the  secondary  user.  The  optimization  problem  in  (3)  is  a 
NP-hard  nonconvex  optimization  problem.  In  the  next  section,  we  propose  a  novel  realizable  suboptimum 
solution  to  this  problem. 


III.  Proposed  scheme 


Using  the  matrix  inversion  lemma,  we  can  express  sfH,  1Sj  as 


/i 


stR_1s-  =  Si  ^  St  i  =  12  k 

’  ’  1-Eih^sfK-W 


A 


Then  (4)  and  (6)  imply  that 


sfR  1Sj  > 


a 


A 


Eihf  +  aEihf 

Thus,  the  optimization  problem  in  (3)  can  be  written  as  follows: 


=  7 i,  i  =  l,2,...,K. 


max  _EssTR,  1  s, 

E„  ss  S  /s 


subject  to  sfH  1s^  >  i  =  1, 2, . . . ,  K 


Ss  Ss  —  1,  Es  ^  En 


(6) 


(7) 


(8) 


3 


Using,  once  again,  the  matrix  inversion  lemma  on  R  1,  i.e. 


TP  ^2t>—  1  c  cT-p—  1 
-p-l  _  T3-1  JZ,sif'ssx'p+nPs^s  -"^p+n 

“  p+n~  i  +  E,Hsm-i„Si  ’ 


(9) 


7s't,sC5s  A^p+n  s 

and  then  combining  (7)  and  (9),  we  can  express  the  optimization  constraints  in  (3)  as  explicit  functions  of 
the  code  sequence  of  the  secondary  user  s.s,  i.e. 


>  EshsSiRp+nssss  Rp+nSj  -in  t,' 

s *  RP+«S*  -  1  +  E  /,2str-  1  s  f  7i’  1  -  1’ 2’  •  •  •  ’  K- 

In  (9),  Rp+n  is  the  autocorrelation  matrix  of  the  primary  signals  and  noise  at  the  base  station,  i.e. 


(10) 


K 


K 


2=  1 


Rp+n  —  yj EihiSibi  +  np)(y^  yj EihiSibi  +  rip)  } 

=  ^2  Eihjsisf  +  cr2 1. 


2=1 


For  notational  simplicity,  let  us  define  the  L  x  L  matrix 


.  —1 


B i  =  hsRp+nstst  Rp+n  -  PihsRp+n 


where  Pi  =  sfRp+ns i  —  -jt.  Then,  the  optimization  problem  (8)  can  be  written  as 


rji  _ -| 

max  x  R7„  x 

X  /s 


subject  to  x1  B?;x  —  Pi  <  0,  i  =  1, 2, . . . ,  K 


(11) 


x'  X  <  E 


—  ^ max 


where  x  is  the  transmitted  signal  vector  of  the  secondary  user,  i.e.  x  =  \/Ejss.  We  note  that  B,.  i  = 
1,2 is  not  positive  semidefinite,  and  the  problem  in  (11)  is  a  nonconvex  quadratically  constrained 
quadratic  program  (nonconvex  QCQP). 

We  observe  that  if  we  use  the  property  7r{AB}  =  7r(BA),  we  are  able  to  represent  the  objective  function 
in  (11)  as 

xTR-1x  =  rr{R-1X},  (12) 
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where  X  =  xx7  .  Thus,  the  optimization  problem  in  (11)  can  take  the  following  form: 

max  rrjR/^X}  (13) 

subject  to  7r{B,;X}  <  fy,  i  =  1,2, ,  K 

7r{X}  <  Emax,  x  y  0,  ranfc(X)  =  1, 

where  X  y  0  denotes  that  the  matrix  X  is  positive  semidefinite.  To  solve  this  problem,  we  propose  to  relax 
the  rank  constraint  in  (13)  and  proceed  by  solving  the  following  problem  instead 

max  7>{R^s'X}  (14) 

subject  to  7r{B,;X}  <  fjl7  i  =  1,2,...,  K 
7r{X}  <  Ernax ,  X  ^  0. 

Then,  (14)  can  be  solved  using  semidefinite  programming.  Strictly  speaking,  we  can  solve  (14)  in  polynomial 
time  only  within  an  error  e  from  the  optimum  solution.  More  specifically,  let  fQ  =  7>-{RT1X}|x=x  where 
X0  is  the  optimum  point,  i.e.  fQ  is  the  optimal  value  of  the  objective  function  in  (14).  Then  for  any  given 
e  >  0,  semidefinite  programming  guarantees  that  we  can  converge  in  polynomial  time  (polynomial  in  the  input 
size  L  and  in  log  1/e)  to  a  solution  that  lies  in  (/<,  —  e,  f0).  In  this  paper,  for  the  semidefinite  programming 
problem  in  (14),  we  propose  to  use  a  primal-dual  interior-point  method  [12],  [13].  In  particular,  we  consider 
the  problem  in  (14)  as  the  primal  optimization  problem,  we  create  a  differently  parameterized  equivalent 
dual  problem  and  then  solve  both  problems  iteratively  in  a  coupled  fashion.  Then,  each  iterations  can  be 
implemented  in  0(L3)  and  the  algorithm  converges  after  a  small  number  of  iteration  to  the  matrix  X*  that 
makes  the  objective  function  Tr{RT1X}  attain  a  value  within  (f0  —  e,  f0).  The  proposed  method  is  outlined  in 
the  appendix.  We  note  that  relaxing  the  rank  constraint  of  the  NP-hard  problem  in  (13)  leads  to  the  optimization 
problem  in  (14)  that  can  be  solved  in  polynomial  time  (by  semidefinite  programming  methods)  as  described 
in  the  appendix.  However,  because  of  the  constraint  relaxation  itself,  the  objective  function  evaluated  at  the 
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convergence  point  of  the  proposed  method  described  above,  7>{R,  1X},  is  just  a  lower  bound  on  the  optimal 
value  of  the  objective  function  in  (13).  To  compute  a  good  feasible  solution  for  our  problem  in  (11),  we 
propose  to  proceed  with  following  randomization  procedure. 

Let  X*  be  the  convergence  point  of  our  proposed  method  described  in  the  appendix.  If  we  select  x  as  a 
Gaussian  random  vector  with  0  mean  and  correlation  matrix  X*,  i.e.  x  /'v(x)  =  A/"(0,  X*),  then  x  is  the 
optimum  solution  of  the  nonconvex  QCQP  in  (11)  “on  average”  over  all  possible  distribution  of  x,  i.e. 

/jv(x)  =  arg  max  ^{x^RT^x}  (15) 

/(*)  ' 

subject  to  E’jx7  B,x}  <  Pi,  i  =  1, 2, . . . ,  K 
T/{x  x}  <  Emax, 

where  /(x)  denotes  the  probability  density  function  of  x.  Implementation  wise,  a  “good”  feasible  vector 
can  be  obtained  by  sampling  x  a  sufficient  number  of  times  and  then,  among  the  feasible  solutions  (i.e.  the 
ones  that  satisfy  the  constraints  in  (15))  we  simply  choose  the  vector  that  maximizes  the  objective  function 

rri  _ 1  i  f  f  m 

X  R/s  x.  However,  global  optimality  can  be  guaranteed  only  if  x  is  sampled  an  infinite  number  of  times, 
which  is  not  realistic.  Our  extensive  studies  indicate  that  if  we  sample  x  until  we  obtain  the  first  feasible 
point,  then  we  can  use  this  vector  as  the  initial  point  of  an  iterative  procedure  that  converges  to  a  good  feasible 
solution.  The  iterative  procedure  is  outlined  below  and  its  performance  is  evaluated  by  simulations  in  the  next 
section.  First  we  express  R/s  as 

R/s  =  SE  ST  +  cr2I,  (16) 

where  S  =  [si, S2, . . . , s^]  denotes  the  matrix  that  has  the  signatures  of  primary  users  as  columns,  and 
S  =  diag(E\q\,  ,  ■  •  • ,  Exq2K).  Using  the  matrix  inversion  lemma,  R  s  can  be  expanded  as 

RT1  =  41  -  4S(S_1  +  STS)”1ST.  (17) 

/s 
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Substitution  of  (17)  into  (11)  leads  to  maximization  of  the  following  objective  function 

“2  x7  x  —  — -j-x^Q  x,  (18) 

aA  cr4 

where  the  matrix  Q  =  S(E'  1  +  SiS)_1Si.  In  (18),  the  first  term  ^x]x  is  a  convex  function  while  the 
second  term  -^xTQx  is  a  concave  function  (the  latter  implies  that  4xTQx  is  convex).  Based  on  the 
first-order  conditions  of  convex  functions  [14],  we  have 

xtx>2x(°)tx-x(0)tx(0),  (19) 

where  x^  denotes  an  initial  feasible  vector.  Then  we  combine  (18)  and  (19)  and  form  an  optimization  problem 

that  maximizes  the  following  concave  function 

Ax(0)Tx  _  1  XTQ  x  _  1  X(0)TX(0)  (20) 

cj4 

that  leads  to  a  suboptimum  solution  for  our  original  problem  in  (11).  To  solve  (20)  we  restrict  all  nonconvex 
constraints  into  convex  sets  (linearization).  In  particular,  we  consider  the  nonconvex  constraints 

xT  B/x  -  fa  <  0,  i  £  lnc,  (21) 

where  lnc  denotes  the  set  of  all  indices  where  x  7  B,x  is  a  nonconvex  function.  Then  we  decompose  the 
matrix  B,  into  its  positive  and  negative  parts  as 

Bt  =  B+  -  B-  (22) 

where  B  =  h‘gTVp+ns,sJ and  B  =  (3ih2s Rp+n  are  all  positive  semidefinite.  Therefore,  the  original 
constraints  (21)  can  be  written  as 

xT B+  x  -  Pi  <  xT B_7  x,  i  €  lnc  (23) 

where  both  sides  of  the  inequality  are  convex  quadratic  functions.  Linearization  of  the  right-hand  side  of  (23) 
around  the  vector  x(0)  leads  to 

xtB+  x  -  Pi  <  x(°)TBr  x(°)  +  2x(°)rBr  (x  -  X(°)),  *  G  Xnc.  (24) 
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In  (24),  the  right-hand  side  is  an  affine  lower  bound  on  the  original  function  xv  B  x.  It  is  thus  implied  that 
the  resulting  constraints  are  convex  and  more  conservative  than  the  original  ones,  hence  the  feasible  set  of  the 
linearized  problem  is  a  convex  subset  of  the  original  feasible  set.  Thus,  by  linearizing  the  concave  parts  of 
all  constraints,  we  obtain  a  set  of  convex  constraints  that  arc  tighter  than  the  original  nonconvex  ones.  Now, 
the  original  optimization  problem  takes  the  form 


max  -lx(°)Tx-  ixTQx-  ix(0)rx(0) 

x  <JZ  <74  GZ 

subjec  to  x3  B+  x  —  x(0)I  B?~  (2x  —  x(0^)  -  /%  <  0,  *  €  Tnc 
X  B;X  f3i  A  0,  i  £  T. nc 
x3x  <  Emax. 


(25) 


The  problem  in  (25)  is  a  convex  QCQP  problem  and  can  be  solved  efficiently  by  standard  convex  system 
solvers  [15]  to  produce  a  new  feasible  vector  x(  1  *  (the  objective  function  in  (11)  evaluated  at  x:  1  j  takes  a 
value  that  is  larger  than  the  value  that  is  based  on  x(0)  y  Repeating  the  linearization  method,  we  can  obtain 
a  sequence  of  feasible  vectors  with  non-decreasing  values  of  the  objective  function  in  (11).  This  procedure 
converges  after  few  iterations.  Our  proposed  scheme  for  power  and  code  allocation  for  the  secondary  user  is 
outlined  in  Fig.  2. 


V.  SIMULATION  STUDIES 

We  consider  a  primary  DS-CDMA  system  with  K  synchronous  users  and  a  pair  of  secondary  transmitter 
and  receiver  with  CDM  mode,  and  the  system  processing  gain  is  L  =  16.  At  the  base  station,  the  transmitted 
SNRs  of  K  users  are  all  equal  to  10 dB  while  the  maximum  transmitted  SNRs  for  the  secondary  user  is  set  to 
8 dB.  All  signatures  for  primary  users  arc  generated  from  the  optimum  binary  signature  set,  which  achieves 
the  Karystinos-Pados  bound  with  the  length  L.  The  SINR  threshold  for  primary  users  is  set  to  3 dB.  The 
channel  coefficients  hn  and  qn,  n  =  1,2 , ,K,  are  considered  as  complex  Gaussian  random  variables  with 


mean  0  and  variance  1,  while  hs  and  qs  are  set  to  0.9.  In  the  randomized  procedure,  the  Gaussian  variable 
is  sampled  1000  times,  and  the  first  and  best  feasible  vectors  arc  denoted  by  x!(,i  and  xmax,  respectively.  In 
the  linearization  procedure,  we  pick  x(°)  as  the  initial  vector,  and  the  output  convergence  vector  is  denoted 
by  xout. 

In  Fig.  3,  we  plot  secondary  transmission  percentage  of  the  SDP  solution  X*  with  rank  1  and  rank  more  than 
1  as  a  function  of  the  number  of  primary  users.  We  observe  that  the  total  secondary  transmission  percentage 
decrease  as  the  number  of  primary  users  increases.  Most  secondary  transmissions  can  be  realized  with  the 
SDP  solution  X*  with  rank  1,  which  means  that  the  optimal  vector  x  can  be  directly  generated  from  the  SDP 
solution  X*  without  the  randomized  and  linearization  procedure.  The  secondary  transmissions  achieved  by 
the  linearization  procedure  arc  below  5%,  especially  almost  zero  when  the  number  of  primary  users  is  greater 
than  19.  It  means  that  most  secondary  transmissions  can  be  realized  without  the  linearization  procedure,  then 
the  large  complexity  of  linearization  is  avoided. 

For  comparison  purposes,  we  evaluate  the  SINR  loss,  SINR(xr_opt )  —  SINR(x),  of  x  =  x^°\  xmax 
and  xout,  with  respect  to  the  SINR  of  the  optimal  real  vector  xr_opt  =  V^max qi,  where  qi  denotes  the 
eigenvector  of  R7S'  with  the  largest  eigenvalue. 

In  Fig.  4,  we  plot  the  SINR  loss  of  x(0\  xmax  and  xoiij  as  a  function  of  the  number  of  primary  users.  Under 
a  limited  number  of  samples,  the  linearized  output  vector  xout  has  a  less  SINR  loss  than  the  “best”  vector 
xmax.  We  observe  that  the  SINR  losses  of  both  xmax  and  xout  arc  much  better  than  x':0:'.  The  linearization 
procedure  improved  the  SINR  performance  from  the  initial  point,  and  meanwhile  realized  the  derandomization. 

V.  CONCLUSION 

In  this  paper,  we  developed  a  power  and  code  allocation  for  a  secondary  user  that  coexists  with  a  primary 
CDM  network.  First,  we  formulated  the  task  as  a  NP-hard  nonconvex  constraint  optimization  problem  (3). 
Then  relaxation  of  the  rank-1  constraint  led  to  a  problem  that  can  be  solved  by  semidefinite  programming. 
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To  do  that,  we  proposed  a  primal-dual  interior-point  method  that  leads  to  a  matrix  solution  (rank  >  1).  A 
subsequent  randomization  procedure  made  it  possible  to  obtain  a  “good”,  on  average,  feasible  solution  after 
searching  over  several  feasible  sampled  solutions.  To  reduce  the  computational  complexity  associated  with  the 
randomization  procedure,  we  developed  an  alternative  method  that  iteratively  solves  the  original  nonconvex 
optimization  problem  by  restricting  all  nonconvex  constraints  into  convex  sets  and  by  initializing  itself  at 
any  feasible  solution  of  the  randomized  procedure  (e.g.  the  first  feasible  sampled  solution  encountered  in  the 
randomization  process). 


APPENDIX 

Interior  point  algorithm 

(i)  Formulate  the  pair  of  primal  and  dual  SDP  problems: 


Primal 

Dual 

max 

X 

THR^X} 

max 
y  z 

bTy 

K 

subject  to 

7>{BjX}  <  Pi,  i  =  1,2, ...  ,K 

subject  to 

y?:Bt  +  yjr+il  =  R^1  +  Z 

i= 1 

7>{X}  <  Ernax .  X  A  0 

y  >  0,  Z^O. 

(ii)  Choose  0  <  <5  <  1  and  define  p  = 

(iii)  Determine  AX,  Ay  and  AZ  by  the  method  in  [13]. 

(iv)  Using  XZ  method,  replace  AX  by  \ (AX  +  AXT). 

(v)  Choose  steplengths  a,  f3  and  update  the  iterates  by 

X  X  +  aAX 

y  <-  y  +  /3Ay 
Z  <—  Z  +  aAZ 
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(vi)  Repeat  (ii)-(v)  until  X  is  feasible  for  the  primal,  (y,  Z)  is  feasible  for  the  dual,  and  the  primal  and 

dual  objective  values  agree  to  a  specified  number. 
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Fig.  2.  The  flow  chart  of  the  proposed  power  and  code  allocation  for  the  secondary  user 
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Secondary  Transmission  percentage(%) 


Fig.  3.  Secondary  transmission  percentage  of  the  SDP  solution  X*  with  rank  1  and  rank  more  than  1  versus  of  the  number  of 
primary  users. 
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Fig.  4.  SINR  loss  of  x-°\  x.max  and  xout  versus  the  number  of  primary  users. 
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